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The electronic structural properties in the presence of constrained magnetization and a charged 
background are studied for a monolayer of FeSe in non-magnetic, checkerboard-, and striped- 
ant if err omagnetic (AFM) spin configurations. First principles techniques based on the pseudopoten- 
tial density functional approach and the local spin density approximation are utilized. Our findings 
show that the experimentally observed shape of the Fermi surface is best described by the checker- 
board AFM spin pattern. To explore the underlying pairing mechanism, we study the evolution of 
the non-magnetic to the AFM-ordered structures under constrained magnetization. We estimate 
the strength of electronic coupling to magnetic excitations involving an increase in local moment 
and, separately, a partial moment transfer from one Fe atom to another. We also show that the 
charge doping in the FeSe can lead to an increase in the density of states at the Fermi level and 
possibly produce higher superconducting transition temperatures. 



Recent experimental advances in molecular beam epi- 
taxy and scanning tunneling microscopy have made it 
possible to study superconducting monolayer systems, 
such as FeSe, which is the simplest iron-based supercon- 
ductor. Studies of FeSe monolayer systems on different 
substrates show significant sensitivity to interface effects 
and give signs of the presence of superconductivity above 
77 K [H, Q . The latter fact is especially interesting be- 
cause bulk samples only show superconducting transition 
temperatures Tc of about 8 K, or 37 K with the appli- 
cation of pressure [sl, 0] • When double-layer graphene is 
used as a substrate, the superconducting gap seems to 
increase with the FeSe film thickness, and a monolayer- 
thick film was found to be non-superconductive [2]. At 
the same time, when FeSe films are studied on SrTiOs, 
scanning tunneling microscopy reveals the presence of 
a 20 meV gap for single unit cell thick films, whereas 
thicker films show no presence of superconductivity [1]. 

Another difference between the bulk and monolayer 
FeSe is the Fermi surface topology: ARPES measure- 
ments show distinctive shape where only pockets at the 
Brillouin zone (BZ) boundaries are present, and no pock- 
ets are seen at the center point F [sj. The same study 
shows the presence of large and nearly isotropic supercon- 
ducting gaps without nodes. Another very recent report 
by the same research group [6] describes the presence of 
two phases with different Fermi surface topology, where 
a superconducting phase can be obtained from a non- 
superconducting one through annealing. Both reports 
give evidences for substantial interface-induced changes 
in electronic structure. Suppression of superconductiv- 
ity by twin boundaries interconnected with the Se-atom 
height with respect to the Fe layer was also found [7]. 

Several theoretical efforts were made to study thin 
films of FeSe [sl, |9|. First-principles study of atomic 
and electronic structures of one- and two- monolayer thick 



films on SrTiOs found semiconductor-like behavior and 
collinear antiferromagnetic order to be present [8[. No 
strong hybridization of the electronic states between the 
substrate and the film was found. Effects of the inter- 
action between the electronic system and soft phonons 
of SrTiOs on FeSe/SrTiOs interfaces was studied using 
functional renormalization group calculations [9]. The 
possibility for the soft phonons at the interface to signifi- 
cantly affect the pairing mechanism and increase Tc with 
respect to the bulk sample was proposed. 

In the current work we study the electronic structure 
of an isolated FeSe monolayer in the non-spin-resolved 
(NM), and for the checkerboard (CH) and the striped 
(STR) antiferromagnetic (AFM) spin configurations. We 
utilize a first principle pseudopotential density functional 
theory and local spin density approximation (LSDA) 
based approach. Our results show that the shape of 
the Fermi surface for the doped CH AFM spin config- 
uration is consistent with the experimentally observed 
one. We study the evolution of the electronic structure 
from the NM to the AFM cases under magnetization con- 
straints on Fe-atoms. Constrained magnetization-based 
estimates show that the transfer of magnetic moment be- 
tween Fe atoms induces a larger bandstructure shift than 
the case when magnetic moments are increased. The 
strength of the magnon-induced energy shifts is up to 
one order of magnitude larger then a typical value for an 
energy shift caused by phonons. We also study the ef- 
fect of the introduction of a charge background, and find 
that it can produce favorable conditions for a higher Tc 
by affecting the number of electrons available for super- 
conducting pairing at the Fermi level. 

The electronic properties are calculated using the gen- 
eralized gradient approximation (GGA) to density func- 
tional theory (DFT)^, [ll| within a planewave pseu- 
dopotential scheme |12l4l4| . We use a 40 Ry cutoff value 



for the kinetic energy of the planewave basis and 600 Ry 
cutoff - for the electronic density. Crystal structures are 
relaxed, so that the forces on the atoms are less then 0.5 
mRy/A. For the NM and the CH spin-polarized con- 
figuration we use a 4-atom tetragonal unit cell. In the 
STR spin-polarized case, we choose a unit cell with 8 
atoms, and the relaxed structure possesses orthorhombic 
symmetry as the striped spin arrangement introduces a 
lattice distortion. The NM, CH and STR phases respec- 
tively show: 6.97, 7.10, 7.11 in atomic units (a.u) for the 
equilibrium lattice constants; 2.57, 2.70 and 2.72 in a.u. 
for the Se-atom heights; and 0, 2.28 and 2.62 bohr mag- 
netons for the equilibrium magnetic moments. The CH 
and STR phases are 13 and 20 mRy, correspondingly, 
lower in total energy then the NM one, when related to 
the 4-atom cell. The background charge is introduced 
using a uniform jellium-like approximation. No relax- 
ation is done for the charged configurations. The magne- 
tization constraints are introduced within LSDA through 
the energy penalty functional, as it is implemented in 
Quantum-ESPRESSO package [l5|, so that a certain value 
of magnetization on the individual atoms is preferred. 

First, we consider the differences in electronic struc- 
ture between the bulk and monolayer cases. As can be 
seen from Fig. HJa), the NM Fermi surface for FeSe is 
very similar to that of the bulk crystal [l6j. There are 
three hole pockets at the zone center F and two elec- 
tron pockets centered at the corners. Fig. [TJc) shows 
the Fermi surface for the STR AFM case, that is similar 
to the bulk case (for bulk FeSe electronic structure see 
[17], for example). There are two pockets that form the 
Fermi surface: the hole-pocket is centered at F, and the 
electron-pocket forms satellites located nearby F towards 
the Y direction. The BZ for the STR case is a fraction of 
those for NM and CH, so we note that the pocket at F for 
the STR configuration's reduced BZ would produce pock- 
ets at both F and M for the other two configuration's. 
The CH case does exhibit a different Fermi surface in the 
monolayer configuration. Fig. [Hb) shows only one elec- 
tronic pocket at M and small hole-like pockets around F 
towards X. The bandstructure analysis shows that this 
difference comes mainly from the change in the shape of 
a "fiat" band formed close to the Fermi level Ep^ below 
it. In contrast to the bulk case, this band shows a signifi- 
cant reduction in energy width. Analysis of the spatially 
resolved integrated local density of states in the region 
of 200 meV below the Fermi level clearly shows the con- 
centration of charge on iron atoms due to the formation 
of the above-mentioned "fiat" band. 

Secondly, since the experimentally observed value 
of the Fe magnetic moment for the iron-based super- 
conductors is usually smaller then the first-principle 
calculations-based theoretical predictions [18|, ll9| , we in- 
troduce constraints on the value of the atomic magne- 
tization to effectively lower the resulting local moment. 
We then study how the electronic structure evolves from 




FIG. 1. (Color online) (a)-(c) Fermi surfaces of an FeSe 
monolayer for the non- spin-polarized, checkerboard AFM and 
striped AFM cases respectively. The non-zero length in z- 
direction is due to the finite size of the unit cell. Brighter 
colored bands are those crossing the Fermi level. NOTE: the 
striped spin-configuration has a different BZ (rotated by 7r/2 
and scaled by ^J2 in kx and /c^-directions). The stripe di- 
rection is toward the given X point, (d) - the unit cell with 
two Fe and two Se atoms, used in our calculations for non- 
magnetic and checkerboard cases, (e) - the unit cell for the 
striped iron spin arrangements; the neighboring atoms are 
also shown. The Fe atoms are red and grey for spin up and 
down, and Se atoms are shown in green (above and below the 
Fe-plane) . 



the non-magnetic to spin-ordered case. Fig. [2] shows the 
resulting evolution for the CH AFM case. It can be seen 
that when we introduce magnetization, the energy states 
split around M, and as the magnetization increases, a 
"flat" band starts forming as a composition of the bands 
both initially present in non-magnetic states around the 
Fermi level and the lower- lying ones. Similar behavior 
is seen at the F point: a state initially at about -0.5eV 
below the Fermi level rises to form the flat band, and the 
hole-like bands lower their positions in energy until they 
fall below Ep- As a result, a sharp peak in the density 
of states is formed right below the Fermi level. Overall, 
the bandstructure undergoes signiflcant change including 
both band rearrangements and shifts. We point out here 
that a similar- looking electronic structure transformation 
is observed in experiment [6] via annealing. 

Since one of the possible effects of the interface could 
be a charge exchange, it is interesting to analyze how the 
electronic structure is affected by doping. To do this, we 
introduce a uniform jellium charge backgroung. Fig. [3] 
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FIG. 2. (Color online) The evolution of the electronic band- 
structure of an FeSe monolayer under constrained magneti- 
zation for the checkerboard antiferromagnetic spin-polarized 
case. From top to bottom the resulting value of the magnetic 
moment on iron atoms is changing from zero (NM - non- 
magnetic) to 2.28 fiB- The units on y-axis are electron-volts, 
and the Fermi level is at zero. To the right the corresponding 
projected densities of states (in states/eV) are plotted for all 
atoms (black solid), Fe-3d-up states (red solid), Fe-3d-down 
states (green dashed) and Se (blue dot-dashed). 
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FIG. 3. (Color online) The electronic bandstructures of the 
FeSe monolayer with different levels of charge doping in the 
checkerboard antiferromagnetic spin-polarized case (with 2.28 
fiB and no constraints). From top to bottom the value of 
excessive background charge is changing from -0.5e to +0.5e, 
as shown. The units on the y-axis are electron-volts, and 
the Fermi level is at zero. To the right the corresponding 
projected densities of states (in states/eV) are plotted for all 
atoms (black solid), Fe-3d-up states (red solid), Fe-3d-down 
states (green dashed) and Se (blue dot-dashed). Red-colored 
(brighter-colored) bands are the ones crossing the Fermi level. 



shows the results for the CH AFM case (with 2.28 jiB 
and no constraints), where two particularly interesting 
things can be noted: first, the width of the "fiat" band 
formed below Ep is very sensitive to positive doping, and 
hole-like states appear at the Fermi level at various point 
of the BZ; secondly, an electron-like state falls below the 
Fermi level at F point for negative doping, forming a sec- 
ond electron-like pocket of the Fermi surface. The latter 
configuration has a Fermi surface that is similar to what 
was obtained for a bulk intercalated FeSe-compound in 
[ioj . The density of states plots given in Fig. |3] exhibit 
a shape, that favors an increase at Ep for both positive 
and negative doping. In the former case the increase is 
induced by moving the Fermi level into the fiat-band- 
originated peak, in the latter case the increase is due 
to the appearance of the above-mentioned electron-like 
band at F. For the NM and the STR phases our cal- 
culations show a rigid band shift oi Ep for the positive 
charge doping. The hole-like bands at F exhibit higher 
sensitivity in their relative position with respect to the 
Fermi level, shifting upward in energy. For the negative 
charge doping we see new electron-like states appear at 



the Fermi level around F for -0.25e in the non-magnetic 
case, and for -0.75e doping in the striped case. The den- 
sity of states at the Fermi level shows a relative increase 
for both doping directions in both NM and STR phases. 

There appears to be a general belief that the non-spin- 
resolved structure studies using the local-density approx- 
imation provide Fermi surface shapes that are in agree- 
ment with experiment at least for 1111 and 122 com- 
pounds [18., ilSj- Nevertheless, some properties such as 
equilibrium lattice constants, for example, are better de- 



scribed within LSDA [2l|, |22|. On the basis of our find- 



ings we suggest here that a specific AFM spin orien- 
tation might be induced in the FeSe layer by the sub- 
strate. There have been previous studies that reported 
that SrTiOa substrate can induce AFM order at the in- 
terface in thin films 0,[2^ . The existence of only one M- 
point centered Fermi surface pocket was addressed with 
respect to the superconducting pairing in [25] and might 
also mean that the currently widely accepted s± model 
[26I . [27! needs to be adjusted. 
To explore the possible pairing mechanisms, we have 
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studied the bandstructure evolution with magnetization 
constraints, and analyze the bandshifts introduced by 
particular magnetic excitations. This approach is analo- 
gous to a frozen phonon approach to evaluate electron- 
phonon coupling. For example, from Fig. [2] we extract 
the information for the relative position of the bottom 
of the electron-like band at M with respect to Ep as a 
function of the Fe magnetic moment. In the unfolded 
non-magnetic BZ of iron lattice only (similar to used in 
[25} ) the introduction of the CH spin orientation corre- 
sponds to the magnon g^- vector of (7r,7r). Similarly for the 
evolution of the NM to the STR configuration we ana- 
lyze the relative position of the hole-like band maxima 
at r. This excitation corresponds to a (tt, 0) magnon in 
the unfolded BZ. Both (7r,7r) and (tt, 0) magnons involve 
an increase of local Fe moments in a corresponding spin 
arrangements with respect to the non-magnetic configu- 
ration. We also apply a small-magnitude magnetization 
shift from (/io, — Mo) to (/io + ^M, —fio-\-Sfi), correspond- 
ing to an effective transfer of local moment between the 
Fe atoms. Our analysis shows that the energy shifts for 
all three cases do not strongly depend on the value of 
the iron-magnetic moment. The relative strengths of the 
shifts are estimated as follows: 

dEM. .dEr. ^dE a m 

'''''^ ' ~di^ ^""'^^ ' ^ = 1 : 1.5 : 6 (1) 

In addition, we analyze the {dii/da) derivative, where a 
is the in-plane lattice constant, to calculate {dE/da)mag 
= {dE / dii){dii/ da). To accomplish this we obtain the 
relaxed structures for each magnetization value consid- 
ered. The estimates give {dE / dV)(^T^^T^^ ^ 1 eV/A, which 
is of the same order as a characteristic value for a frozen 
phonon, if we consider {dE)ph ~ 15 meV, a ~ 3.75A, 
which are typical for FeSe, and a {da/a)ph ~ 0.005. 
So we can conclude that the local moment transfer- 
involving magnetic excitation are coupled to electronic 
states stronger than the excitations that involve magne- 
tization increase only. The relative coupling value for the 
latter is of the same order as for the phonons. 

As seen in the current work, the increase in the su- 
perconducting gap value may originate from the charge- 
transfer-induced increase in the density of states at the 
Fermi level for the FeSe-layer, and the charge transfer 
itself can originate from the interaction with the SrTiOa 
substrate, or from chemical doping in the case of alkali 
doped FeSe materials. Such an increase would mean that 
there are more electronic states that can take part in the 
scattering process forming Cooper pairs, thus increas- 
ing the gap between the superconducting and normal 
states. Controlling the charge carriers density in thin 
films by gating has long been an actively researched field 
for graphene. The results here suggest that similar ap- 
proaches can be beneficial for studying superconductivity 
in FeSe and possibly in other iron-based superconductors. 



Our results also suggest that the superconducting gap in 
an FeSe monolayer can be enhanced with gating. 

In conclusion, we have studied the electronic proper- 
ties of an FeSe monolayer in different spin configurations. 
Our results show that the checkerboard antiferromag- 
netic spin pattern yields a Fermi surface that resembles 
the experimental results for FeSe on STiOs. This sug- 
gests that a similar antiferromagnetic spin pattern is in- 
duced in the monolayer under experimental conditions. 
A study of the evolution of the non-magnetic electronic 
structure into the checkerboard and the striped antifer- 
romagnetic configurations under constrained magnetiza- 
tion reveals that the magnetic excitations involving local 
moment transfer are coupled to electronic states stronger 
then the ones involving an equal change in magnetization. 
By simulating the substrate-induced interface effects we 
show that the presence of a uniform positive background 
charge in the FeSe-layer can lead to an increase in the 
density of states at the Fermi level producing favorable 
conditions for a higher superconducting transition tem- 
perature. Our findings propose that gated devices similar 
to the ones used now for graphene may be able to reveal 
new interesting results in FeSe and other two-dimensional 
superconductors. 
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